Journal  of  Power  Sources  185  (2008)  338-344 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


pri 

SllbiautS 


Nonlinear  model  predictive  control  of  SOFC  based  on  a  Hammerstein  model 

Hai-Bo  Huoa’*,  Xin-Jian  Zhua,  Wan-Qi  Hub,  Heng-YongTua,  Jian  Lic,  Jie  Yangd 

a  Fuel  Cell  Research  Institute,  Shanghai  Jiao  Tong  University,  Shanghai  200240,  China 
b  Institute  of  Process  Engineering,  Chinese  Academy  of  Sciences,  Beijing  100080,  China 

c  School  of  Materials  Science  and  Engineering,  Huazhong  University  of  Science  and  Technology,  Wuhan  430074,  Hubei  Province,  China 
d  School  of  Mechanical  and  Electronic  Engineering,  China  University  of  Geosciences,  Wuhan  430074,  Hubei  Province,  China 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  15  May  2008 

Received  in  revised  form  23  June  2008 

Accepted  25  June  2008 

Available  online  5  July  2008 


Keywords: 

Solid  oxide  fuel  cell  (SOFC) 

Hammerstein  model 

Radial  basis  function  neural  network 

(RBFNN) 

Autoregressive  with  exogenous  input  (ARX) 
Genetic  algorithm  (GA) 

Model  predictive  control  (MPC) 


To  protect  solid  oxide  fuel  cell  (SOFC)  stack  and  meet  the  voltage  demand  of  DC  type  loads,  two  control 
loops  are  designed  for  controlling  fuel  utilization  and  output  voltage,  respectively.  A  Hammerstein  model 
of  the  SOFC  is  first  presented  for  developing  effective  control  strategies,  in  which  the  nonlinear  static 
part  is  approximated  by  a  radial  basis  function  neural  network  (RBFNN)  and  the  linear  dynamic  part  is 
modeled  by  an  autoregressive  with  exogenous  input  (ARX)  model.  As  we  know,  the  output  voltage  of  the 
SOFC  changes  with  load  variations.  After  a  primary  control  loop  is  designed  to  keep  the  fuel  utilization  as 
a  steady-state  constant,  a  nonlinear  model  predictive  control  (MPC)  based  on  the  Hammerstein  model  is 
developed  to  control  the  output  voltage  of  the  SOFC.  The  performance  of  the  MPC  controller  is  compared 
with  that  of  the  PI  controller  developed  in  [Y.H.  Li,  S.S.  Choi,  S.  Rajakaruna,  An  analysis  of  the  control  and 
operation  of  a  solid  oxide  fuel-cell  power  plant  in  an  isolated  system,  IEEE  Trans.  Energy  Convers.  20  (2) 
(2005)  381-387].  Simulation  results  demonstrate  the  potential  of  the  proposed  Hammerstein  model  for 
application  to  the  control  of  the  SOFC,  while  the  excellence  of  the  nonlinear  MPC  controller  for  voltage 
control  of  the  SOFC  is  proved. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Distributed  generation  (DG)  is  a  promising  technology  that  can 
be  used  to  address  some  of  the  technical  as  well  as  environmen¬ 
tal  concerns  in  power  systems.  As  a  kind  of  high-temperature  fuel 
cell,  solid  oxide  fuel  cell  (SOFC)  presents  an  attractive  option  for  the 
DG  technology  because  it  is  modular,  efficient  and  environmentally 
friendly.  Unlike  other  types  of  fuel  cells,  the  SOFC  is  entirely  solid 
state  with  no  liquid  components.  It  usually  works  at  a  high  temper¬ 
ature,  in  the  range  of  800-1000  °C  to  reach  the  electrolytes  ionic 
conductivity  requirement  [1  ]. 

SOFC  is  a  dynamic  device  which  will  affect  the  dynamic  behav¬ 
ior  of  the  power  system  to  which  it  is  connected.  Analysis  of  such 
a  behavior  requires  an  accurate  dynamic  model.  In  the  last  sev¬ 
eral  decades,  fruitful  results  on  modeling  the  nonlinear  dynamics 
of  the  SOFC  have  been  obtained  [2-5].  However,  most  of  these  mod¬ 
els  emphasized  the  detailed  description  of  cell  internal  processes, 
such  as  component  material  balance,  energy  balance  and  electro¬ 
chemical  kinetics,  etc.  These  models  are  very  useful  to  analyze 
the  transient  characteristics  of  the  SOFC,  but  they  are  too  com¬ 
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plicated  to  be  used  in  controller  design.  For  developing  effective 
control  strategies,  Jurado  et  al.  [6,7]  have  presented  identification 
models  for  a  SOFC.  However,  one  of  the  most  important  cell  per¬ 
formance  variables,  fuel  utilization,  has  not  been  examined  when 
the  authors  explored  the  SOFC  dynamic  response  after  the  distur¬ 
bances.  Furthermore,  they  have  not  considered  the  fuel  processor 
in  their  investigation. 

Nonlinear  input-output  block-oriented  models,  such  as  Ham¬ 
merstein  and  Wiener,  do  not  require  much  fundamental  knowledge 
about  a  system,  and  they  are  relatively  easy  to  be  constructed  using 
process  data.  The  Hammerstein  model  consists  of  a  static  nonlin¬ 
ear  block  followed  in  series  by  a  dynamic  linear  block.  It  is  a  type  of 
commonly  used  nonlinear  models,  and  has  been  successfully  used 
to  model  a  class  of  nonlinear  systems  [8-12].  The  identification  of 
the  Hammerstein  system  involves  estimating  both  the  nonlinear 
and  the  linear  blocks  from  the  input-output  observations.  There 
exist  a  large  number  of  identification  methods  for  the  Hammer¬ 
stein  model,  and  most  assume  the  nonlinearity  is  a  polynomial  of 
finite  and  known  order.  However,  if  the  nonlinearity  is  not  a  polyno¬ 
mial  and  the  input  is  not  Gaussian,  these  algorithms  fail  to  converge 
[13].  To  overcome  the  aforementioned  deficiencies,  Narendra  and 
Parthasarathy  [14]  have  pointed  out  that  a  neural  network  could 
be  used  as  a  nonlinear  operator  in  a  Hammerstein  model.  Com¬ 
pared  with  other  types  of  artificial  neural  network  (ANN),  radial 
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Nomenclature 

E0  ideal  standard  potential  (V) 

E  open-circuit  reversible  potential  (V) 

/FC  stack  current  (A) 

J<r  constant  with  the  value  of  N0 /4 F  (mol  s-1  A-1 ) 

/<i  valve  molar  constants  for  hydrogen,  oxygen  and 

water  (mol  s-1  atm-1 ) 

Nf j  hydrogen  reacted  flow  rate  (mol  s_1 ) 

Nf  natural  gas  input  flow  (mol  s-1 ) 

N0  number  of  cells  in  the  stack 

Nfjj1  hydrogen  input  flow  rate  (mol  s-1 ) 

N^  oxygen  input  flow  rate  (mol  s-1 ) 

N°“t  hydrogen  output  flow  rate  (mol  s-1 ) 

Pi  partial  pressures  of  hydrogen,  oxygen,  and  water 

(atm) 

rH_o  ratio  of  hydrogen  to  oxygen 

r  ohmic  loss  (£2) 

T  stack  operating  temperature  (I<) 

u  fuel  utilization 

us  desired  utilization  in  steady  state 

Vdc  stack  output  voltage  (V) 

Greek  letters 

rf  fuel  processor  response  time  (s) 

T[  response  times  for  the  flow  of  hydrogen,  oxygen  and 

water  (s) 


basis  function  neural  network  (RBFNN)  has  a  number  of  advantages, 
such  as  better  approximation  properties,  simpler  network  struc¬ 
tures  and  faster  learning  algorithms.  To  develop  effective  control 
strategies,  a  Hammerstein  model  of  the  SOFC  is  presented  firstly,  in 
which  the  nonlinear  static  part  is  approximated  by  a  RBFNN  and  the 
linear  dynamic  part  is  modeled  by  an  autoregressive  with  exoge¬ 
nous  input  (ARX)  model.  In  addition,  in  the  Flammerstein  modeling 
the  fuel  processor  is  included  and  the  operating  issue  about  the  fuel 
utilization  is  considered  specifically. 

Model  predictive  control  (MPC)  refers  to  a  class  of  algorithms 
that  compute  a  sequence  of  manipulated  variable  adjustments  in 


order  to  optimize  the  future  behavior  of  a  plant.  Garcia  et  al.  [15] 
have  indicated  that  MPC  was  emerging  as  one  of  the  most  pop¬ 
ular  and  effective  control  techniques  in  process  industries.  It  has 
become  an  industry  standard  mainly  due  to  its  intrinsic  capability 
for  dealing  with  constraints  and  multivariable  system. 

SOFC  produces  DC  electric  power  from  fuel  and  oxidant  via 
an  electrochemical  process.  Any  change  in  the  load  circuit  or  its 
demand  for  power  changes  the  operating  conditions  for  the  SOFC. 
As  a  result,  the  output  voltage  of  the  SOFC  will  have  a  high  fluc¬ 
tuation  in  response  to  the  changes.  Furthermore,  if  these  changes 
cause  the  SOFC  to  work  in  the  overused  or  underused  fuel  condi¬ 
tions,  the  actual  performance  of  the  FC  will  be  affected.  In  order  to 
meet  the  voltage  demand  of  DC  type  loads  and  keep  the  fuel  utiliza¬ 
tion  constant  by  controlling  the  natural  gas  input  flow,  MPC  scheme 
based  on  the  Hammerstein  model  and  a  primary  control  loop  are 
developed  in  this  paper. 

The  rest  of  this  paper  is  organized  as  follows.  In  Section  2,  the 
SOFC  dynamic  model  proposed  in  [5,16,17]  is  briefly  reviewed.  The 
detailed  identification  structure  and  identification  algorithms  of 
the  Hammerstein  model  are  given  in  Section  3.  In  Section  4,  the  MPC 
is  formulated.  Some  simulation  results  are  depicted  to  show  the 
validity  of  the  Hammerstein  model  and  the  practical  applicability 
of  the  proposed  control  strategies  to  control  the  fuel  utilization  and 
the  output  voltage  in  Section  5.  Finally,  conclusions  and  suggestion 
for  future  work  are  presented  in  Section  6. 

2.  Theory  for  the  dynamic  model  of  SOFC 

Based  on  the  work  reported  in  [5,16,17],  the  SOFC  dynamic  model 
is  briefly  reviewed  in  this  section.  The  SOFC  dynamic  model  includ¬ 
ing  the  fuel  processor  adopted  in  this  paper  is  shown  in  Fig.  1  [16]. 

2.1.  The  balance  of  plant  (BOP) 

The  BOP  consists  of  the  natural  gas  fuel  storage,  the  fuel  valve 
controlled  by  its  controller,  and  the  fuel  processor  that  reforms  the 
natural  gas  input  Nf  to  the  hydrogen-rich  fuel  N ™ .  In  [5],  the  authors 
introduced  a  simple  model  of  a  fuel  processor  that  converts  fuels 
such  as  natural  gas  to  hydrogen  and  byproduct  gases.  The  model 
is  a  first-order  transfer  function  with  time  constant  i y.  Hence,  the 
fuel  processor  is  simply  represented  by  this  first-order  model. 


Fig.  1.  SOFC  dynamic  model. 
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Although  CO  can  be  a  fuel  in  the  SOFC,  we  suppose  all  CO  will 
take  part  in  the  CO-shift  reaction  if  the  gas  contains  water  [5].  Thus, 
the  overall  cell  reaction  of  the  SOFC  is: 

H2  +  ^02^H20  (1) 

From  Eq.  (1 ),  it  is  seen  that  the  stoichiometric  ratio  between  hydro¬ 
gen  and  oxygen  is  2  to  1.  Oxygen  excess  is  always  taken  in  to  let 
hydrogen  react  with  oxygen  more  completely.  So,  in  this  paper  the 
flow  ratio  of  hydrogen  to  oxygen  is  kept  at  1.145  [16]. 

2.2.  Solid  oxide  fuel  cell 

The  SOFC  consists  of  hundreds  of  cells  connected  in  series  or  in 
parallel.  By  regulating  the  fuel  valve,  the  amount  of  fuel  into  the 
SOFC  can  be  adjusted,  and  the  output  voltage  of  the  SOFC  can  be 
controlled. 

The  Nernst’s  equation  and  Ohm’s  law  determine  the  average 
voltage  magnitude  of  the  fuel  cell  stack.  Flence,  applying  Nernst’s 
equation  and  Ohm’s  law  (taking  into  account  ohmic  losses),  the 
output  voltage  of  the  SOFC  can  be  modeled  as  follows  [16,17]: 


Vdc  =  E-rI¥  c  (2) 

_  . .  _  N0RT .  Ph2Po25 

E  =  N0E0  + In  — - - (3) 

2F  Ph2o 

where 

pH2(s)  =  -  2 I<rhc)  (4) 

Po2 ( s )  =  -  Kr/FC)  (5) 

1  +  Tq2 S  4 

Ph2o(s)  =  1 /^H2°  2/C,/fc  (6) 

1  +  tH20S 

2.3.  Fuel  utilization 


Fuel  utilization  is  one  of  the  most  important  operating  variables 
affecting  the  performance  of  FC.  It  is  defined  as  [5]: 


u  = 


N{?  -N£ut  NL 
h2  h2  _  h2 


n1?2 


nI?2 


Nq/fc 

2FNF1 


(7) 


where  is  the  hydrogen  reacted  flow  rate.  When  the  stack  is 
operated  at  a  high  fuel  utilization,  both  the  voltage  and  power  den¬ 
sity  decrease.  Furthermore,  if  fuel  utilization  is  too  large,  it  becomes 
impossible  for  the  SOFC  to  sustain  the  voltage  across  the  load.  How- 
ever,  it  is  a  great  waste  under  a  low  fuel  utilization  when  there  is  no 
cycling  of  the  anode  gas  flow.  Therefore,  the  fuel  utilization  should 
be  carefully  selected  to  achieve  the  high  SOFC  performance.  In  this 
analysis,  the  fuel  utilization  is  varied  from  0.7  to  0.9,  according  to 
[16,17]. 

To  protect  the  SOFC  stack  and  avoid  large  deviations  in  the  ter¬ 
minal  voltage  due  to  changes  in  the  stack  current,  we  will  hold  the 
utilization  constant.  According  to  Eq.  (7),  the  operation  of  the  SOFC 
stack  with  a  fuel  input  proportional  to  the  stack  current  results  in  a 
constant  utilization  factor  in  the  steady  state.  Thus,  the  SOFC  stack 
is  operated  with  constant  steady-state  utilization  by  controlling  the 
natural  gas  input  flow  to  the  stack  as  [17]: 


Nf  = 


Nq/fc 

2Fus 


(8) 


where  us  is  the  desired  utilization  in  steady  state.  Furthermore, 
because  the  fuel  processor  is  specially  considered  in  the  dynamic 
modeling  of  the  SOFC,  the  relationship  between  a  small  change  of 


Fig.  2.  Hammerstein  model. 


stack  current  A I¥C  and  a  small  change  of  hydrogen  input  ANg1  fed 
to  the  SOFC  stack  can  be  derived  as  [17]: 


ANin  = 

h2  2Fus(1  +  rfs) 


A/, 


FC 


(9) 


3.  Hammerstein  model 


3.1  Problem  statement 


In  this  section,  we  will  consider  the  problem  of  estimating  a 
model  for  a  Hammerstein  system  based  on  the  input-output  data, 

i.e.,  {Uj}1=1 n  and  {y;}i  =  i n.  The  Hammerstein  model  consists 

of  a  RBFNN  for  identification  of  the  static  nonlinearity,  in  series 
with  an  ARX  model  for  identification  of  the  linear  part.  The  struc¬ 
ture  of  the  Hammerstein  model  adopted  in  this  paper  is  illustrated 
in  Fig.  2,  where  u(k)  and  y(k)  are  the  input  and  the  output  of  the 
Hammerstein  model  at  the  kth  sampling  instant,  respectively,  and 
x(/<)  is  the  output  of  the  RBFNN  which  is  usually  unmeasurable.  The 
output  of  the  Hammerstein  model  is  expressed  as: 

Ua  nb 

y(k )  =  -^a,y(k  -  i)  +  ]TV(fe  -j)  (10) 

1=1  j= 0 

where  a*  (i  =  1, . . .,  na)  and  bj  {j  =  0, . . .,  nb)  are  the  parameters  of  the 
ARX  model,  na  and  nb  are  integers  related  to  the  model  order  and 
the  function. 

The  static  nonlinear  part  in  the  Hammerstein  model  is  repre¬ 
sented  by  using  the  RBFNN  depicted  in  Fig.  3  as: 

M 

X(k)=J2wMum  (ii) 

i=i 


where 


4>i  (u(k))  =  exp 


(  llu(fc)  —  c,||2 \ 

V  2df  ) 


(12) 


is  the  Gaussian  function.  M  is  the  number  of  hidden  node,  q  and  d; 
are  the  centers  and  widths  of  the  ith  RBF  hidden  unit,  respectively. 


x(k) 
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Wj  is  the  connection  weights  from  the  ith  hidden  node  to  the  output, 
H  U  denotes  the  Euclidean  norm. 

3.2.  Identification  of  the  Hammerstein  model 

Identification  of  the  Hammerstein  model  involves  estimating 
the  hidden  centers,  the  radial  basis  function  widths  and  the  con¬ 
nection  weights  of  the  RBFNN,  and  the  orders  and  the  parameters 
of  the  ARX  model. 

The  learning  procedure  of  the  RBFNN  mainly  includes  two  parts, 
i.e.,  the  hidden  centers  and  the  radial  basis  function  widths  are 
determined  firstly,  followed  by  the  adjustment  of  the  connection 
weights.  The  standard  RBFNN  training  method  is  time-consuming, 
since  it  requires  examining  many  different  network  structures 
using  a  trial  and  error  procedure.  In  addition,  the  standard  training 
method  determines  the  hidden  centers  by  a  clustering  approach 
which  usually  results  in  a  large  number  of  selected  centers  [18]. 
Furthermore,  it  is  frequently  observed  that  error  correction  learn¬ 
ing  algorithms,  which  aim  to  provide  the  most  appropriate  weights, 
approach  to  local  optima  due  to  their  gradient-based  attitude. 

On  the  other  hand,  in  most  cases  of  the  ARX  identification,  only 
the  parameters  of  the  ARX  model  need  to  be  identified  and  the 
structure  is  given  beforehand.  However,  like  the  other  modeling 
requirements,  the  model  structure  is  a  vital  and  difficult  issue  to 
address  [19].  In  fact,  it  is  hard  to  guess  a  correct  model  structure 
for  the  ARX  model  if  little  information  of  the  system  is  known.  An 
inappropriate  model  structure  will  degrade  the  performance  of  the 
ARX  identification  process  badly  no  matter  which  algorithm  is  used 
to  estimate  the  parameters.  Therefore,  selecting  a  proper  model 
structure  for  the  ARX  model  is  as  critical  as  estimating  the  model 
parameters. 

Genetic  algorithm  (GA)  is  a  kind  of  global  search  algorithm  based 
on  the  mechanics  of  natural  selection  and  population  genetics.  It 
has  been  successfully  applied  to  a  variety  of  optimization  problems. 
For  instance,  it  has  been  employed  to  identify  the  parameters  of  the 
ARX  and  non-linear  ARX  (NARX)  models,  and  to  select  the  optimal 
model  structures  of  RBFNN  and  recurrent  neural  network  (RNN). 
In  this  study,  GA  is  adopted  to  optimize  the  hidden  centers,  the 
radial  basis  function  widths  and  the  weights  of  the  RBFNN,  and 
the  structure  of  the  ARX  model  which  is  entirely  defined  by  the 
integers  na  and  nb  at  the  same  time.  After  the  ARX  model  structure 
is  determined,  the  parameters  of  the  ARX  model  can  be  estimated 
by  the  least  squares  (LS)  algorithm. 

3.2.1  Optimization  of  the  RBFNN  and  the  ARX  model  structure 
based  GA 

GA  is  an  iterative  stochastic  methodology  which  derives  its 
behavior  from  the  process  of  evolution  in  nature.  It  starts  with  a  ran¬ 
dom  population  of  possible  solutions  (chromosomes).  The  fitness 
of  each  chromosome  is  measured  by  computing  the  corresponding 
value  of  a  fitness  function.  Then  new  generations  are  produced  by 
giving  more  probabilities  of  surviving  to  the  individuals  with  the 
best  fitness  values.  As  the  algorithm  proceeds,  the  members  of  the 
population  are  gradually  improved.  The  parallel  searching  mech¬ 
anism  is  the  main  advantage  of  the  GA,  since  it  is  unlikely  to  get 
trapped  in  local  minima. 

Before  optimizing,  some  GA  parameters  such  as  population  size 
P,  crossover  probability  pc,  mutation  probability  pm,  the  maximum 
numbers  of  iterations,  etc.,  need  to  set.  The  optimizing  process  of 
the  GA  is  given  as  follows: 

(1 )  Encoding:  Encoding  is  the  first  and  an  important  part  of  the  GA 
process.  Binary  coding  is  normally  used  and  has  been  proved  to 
be  optimal  in  certain  cases  [20].  In  this  study,  the  binary  coding 
scheme  is  adopted.  In  encoding,  every  chromosome  includes 


the  hidden  centers,  the  radial  basis  function  widths  and  the 
weights  of  the  RBFNN,  and  the  orders  of  the  ARX  model.  To 
obtain  the  desired  accuracy,  eight  bits  binary  code  is  selected 
to  encode  each  parameter  above. 

(2)  Fitness  function:  The  fitness  function  has  a  great  affect  on  the 
convergence  speed  of  the  GA  process.  To  derive  the  GA  search 
process  towards  the  location  of  the  best  solution,  the  fitness 
function  should  be  able  to  reflect  the  key  properties  of  the 
Hammerstein  model.  If  the  fitness  function  contains  inadequate 
information  about  the  Hammerstein  model,  it  will  not  be  capa¬ 
ble  of  identifying  a  chromosome  with  superior  characteristics 
[21  ].  The  fitness  function  in  this  paper  is  chosen  as  [22]: 

fi  —./worst  ~Ji(i  —  1  >  2,  .  .  .  ,  P)  (13) 

where  Jworst  implies  the  Akaike  Information  Criterion  (AIC)  for 
the  worst  individual,  i.e., 

Jworst  —  rnaxjj  (14) 

i 

Here  ,Ji  denotes  the  AIC  of  the  ith  individual  defined  by: 

Ji  =  N  log(Fj)  +  2(na  +  nb)  (15) 

where 

N 

Ei  =  l^(y(k)  -y{i\k))2  (16) 

k=  1 

where,  y^l\k)  is  the  estimated  output  of  the  ith  individual  at 
time  step  “/<”  of  the  Hammerstein  model,  N  is  the  number  of 
sample  points,  £;  is  the  mean  square  error  between  the  actual 
output  and  estimated  output  of  the  ith  individual. 

(3)  Genetic  operations: 

Selection:  After  the  fitnesses  of  all  the  chromosomes  are 
evaluated,  the  GA  enters  the  selection  phase.  According 
to  the  fitness  values,  better  solution  candidates  (individ¬ 
uals)  are  selected.  The  most  familiar  selection  procedure 
is  fitness-proportional  selection  (roulette  wheel  selection). 
So  the  roulette  wheel  selection  strategy  is  adopted  in  this 
study. 

Crossover:  Crossover  is  a  probabilistic  process  that 
exchanges  information  between  two  parent  chromosomes 
for  generating  two  offspring.  Here  single-point  crossover 
with  a  fixed  crossover  probability  of  pc  is  selected. 
Mutation:  Generally,  over  a  period  of  several  generations, 
the  gene  pool  tends  to  become  more  and  more  homo¬ 
geneous  as  one  gene  begins  to  dominate.  So,  a  mutation 
operator  is  introduced  to  guard  against  premature  conver¬ 
gence.  Mutation  randomly  alters  the  gene  from  0  to  1  or 
from  1  to  0  with  probability  pm.  The  mutation  probabil¬ 
ity  is  set  a  very  small  value,  usually  chosen  in  the  range  of 
0.001-0.01. 

3.2.2.  Parameter  identification  of  the  ARX  model 

After  the  ARX  model  structure  is  determined,  the  LS  algorithm 
is  adopted  to  estimate  the  parameters  of  the  ARX  model,  i.e., 

0=  [HrH]_1HTy  (17) 

~  ^  ^  ^  T 

where  0  =  [a\,  a2, . . . ,  ana,  b0,bi, . . .  ,bnb]  is  the  estimated 
parameters  of  the  ARX  model  and 

H  =  [h(l),  h(2), h(JV)]L(na+ni,+l)  (18) 

y  =  [y(i),y(2),  ....yWlLi 


(19) 
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h(k)  =  [-y(k-n-y(k-  2),..., 

-y(k  -  n0),  x(k ),  x(k  -  1 ), . . . ,  x(k  -  nb)]J „Q+nt+1)xl  (20) 

4.  Model  predictive  control 

MPC  refers  to  a  class  of  control  algorithms  in  which  a  dynamic 
model  of  the  plant  is  used  to  predict  and  optimize  the  future  behav¬ 
ior  of  the  process.  The  basic  control  strategy  of  the  MPC  is  to  select 
a  set  of  future  control  horizons  and  minimize  a  cost  function  based 
on  the  desired  output  trajectory  over  a  prediction  horizon  with  a 
chosen  length. 

4.1.  Output  voltage  prediction 

The  MPC  predicts  the  future  output  voltage  of  the  SOFC  with  the 
Hammerstein  model  and  past  input-output  data  in  this  study. 

Based  on  the  established  Hammerstein  model,  the  predictive 
output  voltage  Vdc(k  + 1 )  for  p  steps  ahead  can  be  obtained  by  suc¬ 
cessive  iterations: 

na  nb 

^dc(k  +  0  —  —  ^  \Vdcik  +  l  —  i)+  ^  bjx(l<  +  l  —  j),  0  <  l  <  p  (21) 
i=i  j= o 

where 

xW  =  f„,exp(-!«|^!!!)  (22, 

1=1  '  1  ' 

4.2.  Reference  trajectory 

The  control  goal  of  the  SOFC  is  that  its  output  voltage  can  reach 
the  set  point  quickly  and  smoothly.  The  output  voltage  of  the  SOFC 
should  follow  a  reference  trajectory: 

Vdcr(k  +  j)  =  oJvdc(k)  +  (1  -  a?)V sp  (23) 

where  Vdcr(k+j)  is  the  reference  trajectory  of  the  output  voltage  at 
time  k+j,  Vdc(k)  is  the  actual  output  voltage  at  time  k,  VSp  is  the  set 
point  of  the  output  voltage,  j  =  1, 2, . . .,  p  are  coincidence  points,  and 
a  is  a  parameter  (0  <  a  <  1 ). 

4.3.  Optimization  criteria 

The  objective  of  the  MPC  controller  is  to  compute  future  con¬ 
trol  moves  which  will  minimize  a  criterion  function  based  on  the 
desired  output  trajectory  over  a  prediction  horizon  Hp.  In  predictive 
controllers,  the  predictive  control  law  can  be  yielded  by  minimizing 
the  cost  function.  Hence,  the  choice  of  the  criterion  function  is  of 
paramount  importance  [23].  To  obtain  the  appropriate  control  law, 
the  following  criterion  function  is  minimized.  The  criterion  func¬ 
tion  is  a  quadratic  function  of  the  errors  between  the  predicted 
output  voltage  and  the  desired  trajectory  (reference)  over  the  pre¬ 
diction  horizon,  and  includes  terms  which  penalize  the  excessive 
changes  in  the  manipulated  variables  (the  natural  gas  input  flow). 

Hp  Hc 

J  =  ^[Vdc(k  +  j)  -  Vicr(k  +  j)]2  +  kJ2mk  +  i)  -  Nf(k  +  i  - 1  )]2 
j= 1  1=1 

(24) 

subject  to: 

Nf,  min  —  —  Nf,  max 

where  Vdc{k  +  j)  is  the  predicted  output  voltage  of  the  SOFC,  Hc  is 
the  control  horizon,  Nf(/<  +  i )  is  the  manipulated  variable  at  time  k  +  i, 
and  k  is  the  weighting  factor. 


The  control  signal  (the  natural  gas  input  flow)  changes  only 
inside  the  control  horizon  and  remains  constant  afterward,  i.e., 

Nf(k  +  j)  =  Nf(k  +  Hc  -  1 ),  j  =  Hc, . . . ,  Hp  (25) 

5.  Results 

5.1.  Hammerstein  model  of  the  SOFC 

To  establish  the  desired  Hammerstein  model,  we  choose  the  nat¬ 
ural  gas  input  flow,  the  oxygen  flow,  the  operating  temperature,  and 
the  stack  current  as  the  model  inputs,  the  voltage  as  the  output. 
For  keeping  the  utilization  constant,  the  natural  gas  input  flow  can 
be  controlled  according  to  the  stack  current  which  is  proportional 
to  the  terminal  load.  A  decrease  in  the  load  current  will  in  turn 
increases  the  output  voltage  of  the  SOFC,  while  the  decrease  in  the 
current  decreases  the  natural  gas  input  flow. 

For  the  purpose  of  identification,  the  white-box  model  described 
in  Section  2  is  excited  with  band-limited  white  noise  around  the 
nominal  value  of  the  natural  gas  input  flow  (Nf),  while  all  the  other 
inputs  are  kept  constant  in  their  nominal  values.  The  nominal  oper¬ 
ating  conditions  of  the  SOFC  are  given  in  Table  1  [5,16].  Changes  in 
the  natural  gas  input  flow  are  produced  every  10  s,  with  a  maxi¬ 
mum  amplitude  of  ±70%  of  the  nominal  value.  As  a  result,  a  set  of 
6000  data  points  is  collected  from  the  simulation  with  a  sampling 
time  of  0.1  s.  The  first  3000  data  are  used  to  estimate  the  Hammer¬ 
stein  model  of  the  SOFC,  while  the  following  3000  data  are  used  for 
validation  purposes. 

Through  a  large  number  of  tests,  a  RBFNN  with  6  hidden  nodes 
in  the  Hammerstein  model  is  chosen  for  its  better  performance.  The 
GA  parameters  are  chosen  as:  population  size  P=40,  the  maximum 
number  of  iterations  =  50,  the  total  length  of  every  chromosome 
is  3  x  6  x  8  +  2  x  8,  pc  =  0.6  and  pm  =  0.008.  When  a  load  disturbance 
causes  the  stack  current  to  have  a  step  change  (from  300  A  to  270  A) 
at  150  s,  the  actual  output  voltage  and  the  predicted  output  volt¬ 
age  of  the  Hammerstein  model  are  represented  in  Fig.  4  by  using 
the  above  identification  algorithm.  Fig.  4  indicates  that  the  pro¬ 
posed  Hammerstein  modeling  method  is  applicable  to  describe 
the  nonlinear  dynamic  behaviors  of  the  SOFC,  and  the  Hammer¬ 
stein  model  of  the  SOFC  presented  in  this  paper  is  accurate  and 
valid. 

5.2.  Fuel  utilization  and  output  voltage  control 

According  to  Eq.  (8),  we  can  regulate  the  natural  gas  input 
flow  to  hold  the  utilization  as  steady-state  constant  in  the  pres¬ 
ence  of  external  disturbances.  To  prove  the  validities  of  the  control 
strategies,  we  choose  the  current  disturbance  as  a  multiple  step 


Table  1 

SOFC  operating  point  data 


Item 

Value 

N0 

384 

T 

1273  K 

^FC.rate 

300A 

Us 

0.8 

E0 

1.18  V 

kh2 

0.843  mol  s”1  atirr1 

I< 02 

2.52  mol  s-1  atm-1 

^h2o 

0.281  mol  s-1  atm-1 

th2 

26.1s 

to2 

2.91s 

th2o 

78.3  s 

5s 

r 

0.126  C2 

rH.O 

1.145 
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Fig.  4.  Output  voltage  of  the  actual  and  identified  Hammerstein  models.  Fig.  6>  voltage  control  of  the  SOFC  using  MPC  controller. 


signal  which  reduces  from  300  A  to  270  A  at  150  s,  and  goes  back 
to  300  A  after  650  s.  In  this  situation,  the  control  effect  of  the  fuel 
utilization  is  depicted  in  Fig.  5.  From  Fig.  5,  we  can  see  the  fuel  uti¬ 
lization  of  the  SOFC  can  be  controlled  as  steady-state  constant  by 
regulating  the  natural  gas  input  flow  according  to  the  stack  cur¬ 
rent. 

Based  on  the  Flammerstein  model,  the  MPC  scheme  can  be 
developed  in  succession.  The  objective  of  the  MPC  study  is  to  main¬ 
tain  the  output  voltage  of  the  SOFC  as  a  desired  value  after  the  fuel 
utilization  is  kept  constant.  The  parameters  for  the  MPC  strategy  are 
chosen  as:  Hp  =  40,  Hc  =  6  and  a  =  0.3.  The  proposed  Flammerstein 
model-based  MPC  algorithm  is  simulated  and  the  result  is  shown 
in  Fig.  6,  which  clearly  demonstrates  the  excellent  performance  of 
the  MPC  scheme. 

For  the  purposes  of  comparison,  the  voltage  control  result  of  the 
SOFC  using  PI  controller  proposed  in  [16]  is  depicted  in  Fig.  7.  Com¬ 
paring  Fig.  6  with  Fig.  7,  one  will  notice  that  the  output  voltage  of 
the  SOFC  can  achieve  the  desired  value  no  matter  which  controller 
above  is  adopted.  In  the  case  of  the  above  current  disturbances,  the 
output  voltage  using  the  MPC  controller  can  achieve  the  desired 
value  in  repose.  Flowever,  it  is  more  fluctuant  when  the  PI  controller 
is  used  to  control  the  output  voltage. 


Fig.  5.  Fuel  utilization  control  of  the  SOFC. 


Fig.  7.  Voltage  control  of  the  SOFC  using  PI  controller. 


To  evaluate  the  performance  of  the  control  systems,  the  integral 
of  time  absolute  error  (ITAE)  is  used: 

1200 

ITAE  =  ^k|Vdc(k)-Vsp(k)|  (26) 

k= 1 

The  resulting  ITAE  using  MPC  is  1.4015e+005,  which  is  smaller  than 
that  using  PI  controller  (1.5762e+005).  This  indicates  that  the  MPC 
scheme  based  on  the  Flammerstein  model  has  better  control  effects 
than  that  of  the  PI  controller  developed  in  [16]. 

6.  Conclusions 

To  develop  valid  control  strategies,  a  Flammerstein  model  of 
the  SOFC  has  been  developed  to  describe  the  nonlinear  dynamic 
behaviors  of  the  output  voltage  and  the  natural  gas  input  flow  accu¬ 
rately.  The  Flammerstein  model  adopted  in  this  paper  consists  of 
the  RBFNN  in  series  with  the  ARX  model.  Simulation  results  have 
illustrated  the  applicability  of  the  proposed  Flammerstein  model  in 
modeling  the  nonlinear  dynamic  properties  of  the  SOFC. 

For  protecting  the  SOFC,  a  primary  control  loop  has  been 
designed  to  hold  the  fuel  utilization  as  steady-state  constant.  Then 
based  on  the  Hammerstein  model,  a  nonlinear  MPC  controller  has 
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been  developed  to  control  the  output  voltage  of  the  SOFC  at  a 
desired  value  by  regulating  the  natural  gas  input  flow.  The  good 
control  effect  of  the  fuel  utilization  and  the  favorable  performance 
of  the  MPC  controller  for  voltage  control  of  the  SOFC  have  been 
testified  in  this  study. 

It  is  noteworthy  that  in  the  Hammerstein  modeling  the  oper¬ 
ating  temperature  of  the  SOFC  has  been  assumed  as  constant 
and  the  details  associated  with  concentration  variations  have  not 
been  reported.  In  future  works,  a  control  oriented  SOFC  model 
considering  temperature  effects  and  the  details  associated  with 
concentration  variations  will  be  developed.  Furthermore,  based  on 
this  model,  some  control  strategy  studies  can  be  pursued. 
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